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Resumen 


En este trabajo se han utilizado diferentes procedimientos para estimar 
propiedades hidráulicas de suelos limosos representativos de un sector de la 
llanura argentina, partiendo de información granulométrica y de parámetros 
físico-químicos disponibles en las cartas de suelos de INTA. Se emplean funciones 
de pedotransferencia (FPT) incluidas en los códigos SoilPar (Acutis y Donatelli, 
2003) y Rosetta (Schaap et al., 1999), junto con ecuaciones de regresión entre 
parámetros hidráulicos y parámetros físicos dadas por un conjunto de autores. 
Para el conjunto de las 136 muestras analizadas, la aplicación de los procedimien- 
tos reveló una importante dispersión de resultados. Esto da la pauta de que los 
predictores se deben emplear con precaución, analizando las compatibilidades 
entre las muestras de suelo de estudio y las que se utilizaron originalmente en 
la propuesta de las ecuaciones. Se estimaron los parámetros medios ponderados 
para la cuenca del A” Ludueña y los pronósticos se compararon con valores 
obtenidos mediante calibración de un modelo matemático hidrológico. A partir 
de los resultados obtenidos pudo concluirse que el código que mejor se aproximó 
a los parámetros hidráulicos calibrados fue el programa Rosetta, suministrándole 
la información completa (composición textural, densidad aparente y dos puntos 
de la curva de retención). Eventualmente pueden utilizarse las formulaciones de 
Saxton et al. (1986) y Wosten et al. (1999) para estimar la conductividad hidráulica 
saturada y el parámetro a. de la ecuación de van Genuchten, respectivamente. Los 
resultados alcanzados al presente son promisorios y validan esta metodología de 
ponderación sencilla para ser aplicada en otras cuencas de la región. 


Palabras clave: suelos limosos, parámetros hidráulicos, funciones de pedotrans- 
ferencia, llanura argentina. 


Introducción tecnología actual de modelos sofisticados, 
que se encuentra en un estado avanzado con 
El interés en la zona no saturada (ZNS) ha respecto al conocimiento de los parámetros 
aumentado en los últimos años debido a la que engloban esos modelos. En trabajos de 
evidencia creciente de la afectación ambiental calibración de modelos hidrológicos continuos 


provocada por actividades agrícolas e industria-— se ha demostrado la alta sensibilidad que pre- 


les. En este ambiente, los modelos numéricos 
son rutinariamente usados en investigación y 
gestión para predecir el movimiento de agua 
y solutos en la ZNS. Los resultados de las 
aplicaciones sólo pueden ser confiables si las 
propiedades del medio poroso son conocidas. 
Esto representa la principal limitación a la 


sentan los procesos subterráneos, observables en 
el sistema mediante los niveles freáticos, y los 
procesos hidrológicos superficiales, observables 
a través de volúmenes de escurrimiento, frente 
a los parámetros hidráulicos que caracterizan 
la zona no saturada (Zimmermann, 2003). Esto 
resalta el rol protagónico que ejerce la zona 


Tecnología y Ciencias del Agua, antes Ingeniería hidráulica en México, vol. IL, núm. 1, enero-marzo de 2011, pp. 99-116 


Ciencias del Agua, vol. Il, núm. 1, enero-marzo de 2011 


« vj30]0u72], 


Zimmermann y Basile, Estimación de parámetros hidráulicos en suelos limosos mediante diferentes funciones de pedotransferencia 


vadosa en la recarga freática y en la infiltración, 
que constituyen los procesos que interconectan 
la hidrología superficial y la subterránea. 
La zona no saturada debe ser analizada con 
detenimiento, ya que en ella se da lugar a 
importantes procesos de intercambio de flujo. 

Para caracterizar el medio no saturado se 
requiere del conocimiento de las curvas de 
conductividad k(0) y retención wy(0) o bien, 
difusividad D(0), siendo 8 el contenido 
volumétrico de humedad, k(0) la conductivi- 
dad hidráulica no saturada, y(0) el potencial 
mátrico del suelo y D(0) el coeficiente de 
difusividad. 

En trabajos previos (Zimmermann 2006a, 
2006b) se han estimado para los suelos típicos 
de la llanura pampeana santafecina valores de 
referencia de los parámetros que gobiernan 
sus funciones hidráulicas. Estos valores fueron 
obtenidos de distintas fuentes de información: 
calibraciones de modelos matemáticos hidro- 
lógicos, resultados de ensayos compilados de 
estaciones del Instituto Nacional de Tecnología 
Agropecuaria (INTA) y las curvas de retención 
calibradas con la información de la base de 
datos pública recopilada por el Departamento 
de Agricultura de Estados Unidos, USDA, 
denominada UNSODA (Leij et al., 1996). En 
el mismo trabajo se propusieron regresiones 
entre parámetros hidráulicos y parámetros 
físicos (densidad húmeda, contenido orgánico, 
porosidad y pH), con el fin de pronosticar los 
primeros en función de los demás parámetros. 

Utilizándose la misma base de datos 
UNSODA se han propuesto regresiones lineales 
y no lineales entre parámetros hidráulicos de 
la curva de retención, y parámetros físicos 
y granulométricos de las muestras de suelo 
(Zimmermann y Basile, 2006, 2007). 

Existen otras funciones de pedotransfe- 
rencia, entre las cuales pueden citarse las 
propuestas por Rawls et al. (1982), Rawls y 
Brakensiek (1985), Cosby et al. (1984), Saxton et 
al. (1986), Vereecken et al. (1989, 1990), Sharpley 
y Williams (1990), Hutson y Wagenet (1992), 
Wósten et al. (1999) y otros. También se han 
propuesto códigos informáticos tales como 


el SoilPar (Acutis y Donatelli, 2003) y Rosetta 
(Schaap et al., 1999), que incluyen relaciones 
similares que han sido ajustadas a otras bases 
de información edafológica. 

En este trabajo, el conjunto de funciones 
de pedotransferencia mencionadas preceden- 
temente se han aplicado a perfiles tipo de suelos 
de la provincia de Santa Fe (área pampeana, 
Argentina) y se han comparado, con el objetivo 
de analizar su semejanza en los resultados. 


Modelos analíticos para curvas de 
retención de humedad y conductividad 
hidráulica no saturada 


Aunque se han propuesto muchas funciones 
empíricas para la curva de retención son pocas 
las que tienen una base fenomenológica. Una 
de las más conocidas y populares es la de 
Brooks y Corey, según la cual: 
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donde 6, es el contenido de humedad para 
mediosaturado; 0, lahumedad correspondiente 
a saturación irreductible; A, un parámetro que 
depende de la distribución de los tamaños de 
poros; Y, la succión crítica (presión de entrada 
de aire), y S, es la saturación efectiva. Para 
soslayar las limitaciones de la curva de Brooks 
y Corey, específicamente la discontinuidad en 
la zona de saturación, se han propuesto curvas 
con forma sigmoidal. Entre ellas se destaca 
la de van Genuchten, cuya expresión es la 
siguiente: 


S.= [1 + (owy)"] 0 (2) 


donde Q, m y n son parámetros empíricos. 
Este modelo de curva de retención tiene 
más flexibilidad para adaptarse a las curvas 
reales, ya que depende de cinco parámetros 
independientes (a, m, n, O, 0). En la práctica 
resulta conveniente restringir los valores de 
m y n, de forma que m = 1 - 1/n. Con ello se 
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obtienen expresiones sencillas de la curva de 
conductividad que se deriva de ella. La curva 
de van Genuchten incluye la de Brooks y 
Corey como un caso particular cuando n tiende 
a infinito, siendo el producto n.m constante. 
En este caso, A es igual a 1.m y a=1/w.. 

Para la estimación de la conductividad 
hidráulica relativa, K_= K (0) / K, siendo K. 
la conductividad hidráulica saturada; con 
el contenido de humedad suelen utilizarse 
fórmulas empíricas del tipo siguiente (Irmay, 
1954): 


Este tipo de curva suele conducir a buenos 
ajustes en suelos de granulometría uniforme. 
Otros autores, como Mualem (1976), proponen 
expresiones de K_ obtenidas a partir de la curva 
de retención. Utilizando la curva de retención 
de van Genuchten se obtiene una expresión 
para K, en el caso particular en el que m = 1 — 
1/n, K se reduce a: 


K,=S." [1- (1-5, my")? (3) 


donde el parámetro p asume el valor de 0.5. 


Parámetros para las funciones 
hidráulicas de suelos no saturados 


Existen métodos indirectos para generar las 
propiedades hidráulicas de suelos, que varían 
en términos de metodología y complejidad. 
Pueden distinguirse tres grupos principales: 
los métodos basados en la distribución del 
tamaño de poros, los métodos inversos y las 
funciones de pedotransferencia. 

Los primeros se usan frecuentemente 
para estimar conductividad hidráulica no 
saturada y la curva de retención de humedad 
con las funciones de Brooks-Corey o Van 
Genuchten. 

Se han propuesto métodos cuasi-físicos 
que usan el concepto de similitud de la forma 


entre distribuciones del tamaño de poros y de 
la partícula (Arya y Paris, 1981; Haverkamp y 
Parlange, 1986). 

Los métodos inversos están basados en 
la resolución numérica de la ecuación de 
Richards junto a un algoritmo de optimización 
y paralelamente mediciones de campo o 
del 
modelo a los medidos se pueden obtener los 


laboratorio. Ajustando los resultados 


parámetros hidráulicos del suelo en forma 
indirecta. Los métodos inversos son a menudo 
vulnerables por la no singularidad de los 
resultados, es decir, dos o más juegos de 
parámetros optimizados pueden ser aplicables 
al problema estudiado. 

Las funciones de pedotransferencia (FPT) 
ofrecen un tercer método para estimar 
propiedades hidráulicas, usando el hecho de 
que las mismas se relacionan con la compo- 
sición textural del suelo entre otra información 
taxonómica (por 
distribución de tamaño de partículas, densidad 


disponible ejemplo, la 
aparente, contenido de materia orgánica, 
etcétera). Sin embargo, la gran mayoría de 
FPTs son completamente empíricas y no usan 
ningún concepto físico. En este trabajo se ha 
empleado esta tercera vía para la estimación de 
los parámetros hidráulicos. 

Funciones de pedotransferencia utilizadas 

Existe un conjunto de ecuaciones que 
relacionan propiedades hidráulicas de suelos 
con parámetros físicos y granulométricos de 
las muestras. Un primer tipo de FPT se basa en 
clasificar los suelos por su composición textural 
y asignar iguales parámetros, asumiendo que 
suelos similares tienen parámetros hidráulicos 
similares (Carsel y Parrish, 1988; Rawls et al., 
1982). 

Otro grupo de FPT se basa en ecuaciones 
de regresión, ofreciendo una solución con 
propiedades hidráulicas que varían en forma 
continua a través del triángulo de clasificación 
textural del USDA. Las predicciones pueden 
ser mejoradas usando propiedades físicas del 
suelo como información adicional, tales como 
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la densidad aparente, porosidad o contenido 
de materia orgánica. Otros autores mejoran 
la predicción incluyendo uno o más puntos 
de la curva de retención (Rawls et al., 1992; 
Williams et al., 1992). Se ha comprobado que 
este tipo de FPTs mejora la predicción de los 
parámetros hidráulicos de los suelos de la 
zona de estudio (Zimmermann y Basile, 2008). 
Tomando en cuenta esta concepción previa, 
se ha ampliado el grupo de FPTs a considerar, 
y en el cuadro 1 se muestran las ecuaciones 
de regresión que se han empleado para las 
comparaciones en este trabajo. En la misma 
se incluyen regresiones lineales múltiples 
entre los parámetros hidráulicos y otras 
características físicas y granulométricas para 
muestras de suelos limosos extractadas de la 
base UNSODA (Zimmermann y Basile, 2006, 
2007). 

Para la implementación práctica de algunas 
FPT se han desarrollado varios códigos de 
cálculo, entre los que pueden mencionarse 
los programas RosettaOd y SoilParO. Ambos 
fueron utilizados en este trabajo y se describen 
sucintamente a continuación. 


Programa Rosetta 


Rosetta implementa funciones de pedotrans- 
ferencia que predicen las curvas de retención 
y conductividad basadas en ajustes de 
composición textural, densidad aparente, y 
en uno o dos puntos de la curva de retención 
(punto de permanente y/o 
capacidad de campo). Rosetta fue desarrollado 
por Marcel G. Schaap en el United States 
Salinity Laboratory (USDA), 
California, Estados Unidos, 

por Leij y van Genuchten. Este programa 


marchitez 


Riverside, 
supervisado 


emplea para el ajuste de parámetros redes 
neuronales para mejorar las predicciones 
de las FPTs empíricas. Los parámetros 
hidráulicos óptimos, relacionados con datos 
de la entrada (tamaño de partícula, densidad 
aparente, etcétera) se obtienen mediante un 
procedimiento de calibración iterativo. El 


software es de distribución libre y gratuita 


información 
consultando la página web http: / /www.ussl. 
ars.usda.gov. 


y puede obtenerse mayor 


Programa SoilPar 


SoilPar 2.0 es un programa para estimar 
parámetros hidráulicos del suelo desarrollado 
por Marco  Acutis (Departamento de 
Ingeniería Agrícola y Agronomía Territorial, 
Nápoles, Italia) y Marcello  Donatelli 
(Instituto de Investigación para Cultivos 
Industriales, Bologna, Italia). El programa 
permite guardar información de suelos en 
un banco de datos georreferenciado, estimar 
parámetros hidráulicos de suelos usando 
varios procedimientos incorporados al 
programa, comparar las estimaciones contra 
datos medidos en forma gráfica y/o tabulada, 
y crear mapas que usan el formato ESRI. 
Once métodos estiman uno o más de 


las características siguientes: contenido 
de humedad para diferentes tensiones 
matriciales, conductibilidad hidráulica 


saturada y densidad aparente. Tres métodos 
estiman los parámetros de curvas de retención 
(Brooks y Corey, 1964; Hutson y Cass, 1987; 
van Genuchten, 1980), y un método para 
estimaciones de conductibilidad hidráu- 
lica no saturada y parámetros de curvas 
de retención (Campbell y Shiozawa, 1994). 
El software corre bajo Windows 98/NT/2000/ 
XP y es de acceso libre vía Internet. Los 
procedimientos de pedotransferencia se 
clasifican en: a) estimación puntual de algunos 
valores específicos de interés de la curva 
de retención de humedad, conductibilidad 
hidráulica no saturada y densidad aparente; 
b) estimación de parámetros de curvas de 
retención. El programa permite el ajuste de 
parámetros para las curvas de retención en 
caso de contar con datos medidos de contenido 
de humedad y potencial matricial mediante 
un procedimiento de calibración no lineal. 
El mismo está basado en el método simplex 
con restricciones, para evitar inconsistencias 
físicas en los valores ajustados. SoilPar es un 
software disponible gratuito con propósitos 
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Cuadro 1. Funciones de pedotransferencia empleadas en este trabajo. 


Ml Funciones de pedotransferencia Fuentes 


8,=0.1 
8, =0.7919 + 0.001691.C — 0.29619.DA — 0,000001491.L?+ 0,0000821.M0*+ 0.02427.C“ 
+ 0.01113.L” + 0.01472.In(L) - 0,0000733.MO.C - 0.000619.DA.C - 0.001183.DA.MO 
-0.0001664.top.L 


14.96 +0.03135.C + 0.0351 .L +0.646.MO +15.29.DA —0.192.top - 4.671.DA? Wósten et al., 1999 
-0.000781.C?- 0.00687.M0?+ 0.0449. MO”? + 0.0663. In(L) + 0.1482 .In(MO) 


-0.04546.DA.L — 0.4852 .DA.MO + 0.00673.top.C 
0.=e€ 


-25.23 -0.02195.C +0.0074.L— 0.194.MO + 45.5.DA —7.24.DA? + 0.0003658 .C? 
+0.002885.MO?- 12 .81.DA7*-0.1524.L?-0.01958.MO”? - 0.2876 .In(L) 
0.0709. In(MO) — 44.6. In(DA) —0.02264.DA.C +0.0896 .DA.MO +0.00718 .top.C 


n=1+e 

0, = 0.015 + 0.005.C + 0.014.CO 

6, = 0.803 — 0.283.DA + 0.0013.C Vereecken et al., 
In(a) = 2.486 + 0.025.A — 0.351.C 1989 

In(1) = 0.035 — 0.009.A — 0.013.C + 0.015.A? 


1 


5.3396738 + 0.1845038.C — 2.48394546.0, — 0.00213853.C2 — 0.04356349.A.0s 
-0.61745089.C.9s + 0.00143598.A?.0% — 0.00855375.C?.8? — 0.00001282.A?.C 
+0.00895359.C2.8; — 0.00072472.A?.0s + 0.0000054.C2.A + 0.5002806.02.C Rawl 
e aw y 


Brakensiek, 1985 
—0.7842831 + 0.0177544.A — 1.062498.0, — 0.00005304.A?2 — 0.00273493.C2 
+1.11134946.0? — 0.03088295.A.0, + 0.00026587A?.0? — 0.00610522.C282 


—0.00000235.A2.C + 0.00798746.C2.8, - 0.00674491.0?.C 
n=1+e 
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Textura arcilla limosa 
6. = 0.8757 + 0.0847n — 0.31880, 
0. = 0.1264 — 0.07778, — 0.00334C0O + 0.11650, 


n = 1.1466 + 1.31610, 

Textura franco arcillo-limosa 
9, =0.1261 + 0.4006n - 1.273, Zimmermann y 
o. = 0.2117 — 0.1444n — 0.05356, + 0.001325pH — 0.0218, Basile, 2007 
n = 0.7971 + 0.9498n 

Textura franco-limosa 

6. = 0.5318 + 0.1578n — 0.1528, + 0.0003682CVg 
0.= 0.1888 — 0.08636, + 0.003459C0O + 0.1170, 
n = 1.56 — 0.11830, — 3.856501 


0, = (50.5 — 0.142.A — 0.037.C)/100 Cosby et al., 1984 


6, = 0.332 — 0.0007251.A + 0.1276.10g,, (C) Saxton et al., 1986 
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Cuadro 1. Funciones de pedotransferencia empleadas en este trabajo. 


o Funciones de pedotransferencia Fuentes 


K¿= 1019.6-081 log(L)- 1.09 log(C) - 4.64.DA] 


Jabro, 1992 


K¿ = 2.94.101-0.6 + 0.0126.A—0.0064.C] Cosby et al., 1984 
K 1 


24.10 (0.044 A — 0.61.DA) 


Jaynes y Tylor, 1984 


Ks = 1 570 e-0.1975.C Puckett et al., 1985 


Kgs = 0.54.e(-0.007.A — 0.167.C) 


s=e 


K;, = 8.29 - 0.0782.C + 0.085.L — 4.73.DA 


K¿=0.04167,el20.62- 0.96 In(C) - 0.66 In(A) - 0.46 In(MO) - 8.43 DA] 
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Ks= e[12.012 — 0.0755.A +[-3.895 + 0.03671.A — 0.1103.C + 0.00087546.C?) / 8] 


7.755 + 0.0352 L + 0.93 top — 0.967 DA? — 0.000484 C2 — 0.000322 L2 + 0.001 L”1 


Campbell y Shizawa, 
1994 


Saxton et al., 1986 


Rawls y Brakensiek, 
1985 


Gulser y Candemir, 
2008 


Vereecken et al., 1990 


Wosten et al., 1999 


— 8.96847 + 19.52348.8s — 0.028212.C + 0.00018107.A2 — 0.0094125.C2 
- 8.395215.0% + 0.077718.A.0s — 0.00298.A2.02 ss 0.019492.C2.0% 
+ 0.0000173.A2C + 0.02733.C2.8, + 0.001434. A2.8, — 0.0000035.C2.A 


0.0748 MO”! - 0.643 In(L) — 0.01398 DA.C — 0.1673 DA.MO + 0.02986 top.C 


[e =0.04167. e —0.03305 top.L 


K, = 0.04167,e/-42.6 + 8.71 MO + 61.9 DA-20.79 DA?- 0.2107 MO? - 0.01622 CMO - 5.382 DA MO] 


K;¿=1 058.5(0, — 0_33)9355 


no comerciales. El paquete de la instalación 
puede obtenerse de http: / /www.isci.it /tools. 


Aplicación a las series de suelos típicas 
en la región sur de la provincia de Santa 
Fe (pampa húmeda argentina) 


La llanura chaco-pampeana de Argentina, de 
clima templado húmedo, está caracterizada 
por un relieve de moderada pendiente 
topográfica. El área meridional de la 
provincia de Santa Fe, motivo de interés en 
este estudio (figura 1), presenta un módulo 
pluviométrico del orden de los 1 000 mm y 
la temperatura media anual del orden de 
los 17 *C. La unidad estratigráfica superior 
del paquete sedimentario, que constituye 


el sustrato geológico de los sistemas 


Wosten et al., 2001 


Ahuja, 1989 


hidrológicos, está compuesta por el loess de 
la formación pampeana, integrada por limos 
y arcillas de colores castaños. El espesor de la 
formación pampeana, en el sector de análisis, 
varía entre 30 y 150 metros. Los mapas de 
suelos disponibles en la provincia de Santa 
Fe, elaborados por el INTA, clasifican, desde 
un punto de vista textural, los suelos de la 
región entre las fracciones franco limosa, 
franco arcillo limosa y arcilla limosa. 

A fin de aplicar los procedimientos descri- 
tos, se contó con información elaborada por 
INTA, que caracteriza los perfiles de suelos 
presentes en el área de estudio. De las cartas de 
suelos disponibles se extractó la información 
de los perfiles tipo siguientes: A. Juárez, 
Arequito, Casilda, Chabas, Faustina, Gelly, 
Hansen, Huemul, Juncal, Las Gamas, Lima, Los 
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Figura 1. Área de estudio. 


Nogales, Maizales, Monte Fiore, Pergamino, 
Peyrano, Roldán, Santa Lucía, Villada, Villa 
Eloisa y Zavalla. 

En cada uno de los 21 perfiles, en función de 
su profundidad y su estratigrafía, se detallan 
entre tres y ocho subhorizontes, de los cuales las 
cartas de suelo publicó la granulometría, junto 
a otros parámetros biofísicos, como contenido 
de materia orgánica, capacidad de intercambio 
catiónico, pH, etcétera. El conjunto de estratos 


para todas las series estudiadas totalizan 136 
muestras de suelo. A manera de ejemplo, 
en el cuadro 2 se presentan los datos y las 
estimaciones realizadas con el software SoilPar 
para la serie Casilda, con sus ocho estratos. 

El mismo trabajo se realizó para cada uno 
de los horizontes de las series, obteniéndose 
densidades aparentes, humedades para capa- 
cidad de campo (CC), punto de marchitez 
permanente (PMP) y otras tensiones matri- 
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Cuadro 2. Propiedades granulométricas, físicas e hidráulicas estimadas para la serie Casilda. 


Unidades L1 1É2 L3 L4 L5 L6 L7 Ls 

Profundidad (m) 0.18 0.27 0.48 0.76 1.04 1.28 1.80 2.10 
Densidad seca (Rawls) (t/m3) 1.46 1.44 1.42 1.38 ES 9 1.42 1.45 1.44 
[PMP (Brakensiek y Rawls) (m/m) 0.19 0.22 0.22 0.28 0.23 0.19 0.14 0.16 
[CC (Brakensiek y Rawls) (m/m) 0.39 0.42 0.40 0.46 0.40 0.35 0.31 0.32 
Contenido de arena muy gruesa (%) 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 
¡Contenido de arena gruesa (%) 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 
¡Contenido de arena media (%) 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 
Arena fina (%) 2.00 1.50 1.00 1.00 1.50 1.00 1.80 1.30 
[Arena muy fina (%) 8.20 5.50 4.20 2.50 4.60 6.90 9.80 6.90 
¡Contenido de limo (%) 69.90 64.40 61.70 49.60 55.80 62.80 67.80 68.20 
¡Contenido de arcilla (%) 19.87 28.57 33.07 46.87 38.07 29.27 20.57 23.57 
Coarse fragment content (%) 

Materia orgánica (%) 2.55 2.29 1.22 0.99 0.49 0.27 0.15 0.10 
Capac. intercambio catiónico (cmol /kg) 16.70 19.30 21.40 30.60 29.80 24.70 18.80 20.60 
[pH - en agua 6.20 6.40 6.60 7.00 7.10 6.80 7.50 8.20 
Humedad a 10 (-J/kg) (m/m) 

¡Humedad a 20 (-J/kg) (Brakensiek y Rawls) (m/m) 0.43 0.46 0.44 0.49 0.44 0.39 0.35 0.36 
Humedad a 30 (-J/kg) (Brakensiek y Rawls) (m/m) 0.39 0.42 0.41 0.46 0.40 0.35 0.31 0.32 
[Humedad a 50 (-J/kg) (Brakensiek y Rawls) (m/m) 0.34 0.38 087% 0.42 0.36 0.32 0.27 0.28 
¡Humedad a 100 (-J /kg) (Brakensiek y Rawls) (m/m) 0.31 0.34 0.33 0.39 0.33 0.29 0.24 0.26 
[Humedad a 500 (-J/kg) (Brakensiek y Rawls) (m/m) 

¡Humedad a 1 500 (-J /kg) (Brakensiek y Rawls) (m/m) 0.19 0.22 0.23 0.28 0.23 0.19 0.15 0.16 
¡Conductividad hidráulica saturada (Puckett) (mm/h) 31.01 5.56 2.29 0.15 0.85 4.84 27.00 14.93 
van Genuchten al (1/cm) 0.0014 0.0010 0.0017 0.0010 0.0019 0.0036 0.0051 0.0053 
ivan Genuchten m (adim) 0.53 0.19 0.17 0.14 0.41 0.35 0.37 0.33 
van Genuchten n (adim) 2.14 1.24 IL 1.16 1.68 1.54 1.58 1.48 
van Genuchten 0, (adim) 0.18 0.01 0.03 0.02 0.20 0.16 0.13 0.13 


Referencias: los valores en itálica son estimados aplicando el programa SoilPar. 


ciales, conductividad hidráulica saturada y 


parámetros de van Genuchten, para la curva 


de retención de humedad. 


También se aplicaron las ecuaciones de las 


regresiones multivariadas presentadas en el 


cuadro 1. 


Se usó el software Rosetta bajo tres diferentes 


grados de información de entrada: (a) com- 
posición textural (arena, limo y arcilla, SSC), 
(b) composición textural y densidad aparente 
(SSCB) y (c) lo anterior, junto a humedades 
de CC y PMP estimadas. Para la clasificación 
textural de las muestras se consideraron los 
contenidos de arena, limo y arcilla disponibles 
como información en la base del Instituto 
Nacional de Tecnología Agropecuaria. 


Resultados y discusión 


Los resultados de la aplicación de estos criterios 
conjuntos se muestran en las figuras 2 a 6. 
Como se aprecia en las figuras, las diferentes 
metodologías empleadas arrojan resultados 
disímiles, especialmente en la estimación de 
la conductividad hidráulica saturada K. En 
trabajos previos (Zimmermann y Basile, 2008) se 
compararon los resultados de SoilPar y Rosetta 
para un conjunto de 321 muestras representa- 
tivas de suelos en la región de estudio, los que 
fueron clasificadas por texturas. Se concluyó 
que los resultados obtenidos por Rosetta 
fueron los que mejor se acercaron a valores 
calibrados mediante modelación matemática 
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hidrológica realizada en la cuenca del arroyo 
Ludueña, ubicada en la región de estudio. Por 
ello, en las comparaciones realizadas en este 
trabajo, se han tomado como referencia los 
parámetros estimados con el programa Rosetta. 
Allí se atribuyeron las texturas franco limosas 
para caracterizar el horizonte “A” y las arcillo 
limosas para el horizonte “B”. 

La figura 2 muestra las estimaciones 
de la humedad residual. Allí se observa 
que cualquiera de las opciones de Rosetta 
pronostica 6, dentro de un rango de variación 
aproximado al 20%. El programa SoilPar 
presenta patrones disímiles y evalúa tanto por 
defecto como por exceso. La formulación de 
Vereecken (1989) propone estimaciones por 
exceso e incluso valores bastante alejados de 


los promedios recomendados por diferentes 
fuentes bibliográficas. 

La figura 3 muestra las estimaciones 
para la humedad de saturación. En términos 
generales, hay mucho menos dispersión 
que para la humedad residual. Casi todos 
los pronósticos caen dentro del +20%. Se 
podría afirmar que es el parámetro que 
puede pronosticarse con menor variabilidad. 
Debido a que la información granulométrica 
disponible en las cartas de suelo no incluye una 
curva de distribución, no fue posible aplicar 
las ecuaciones propuestas por Zimmermann y 
Basile (2007). 

En la figura 4 se observa que las tendencias 
en el pronóstico del parámetro A. son muy 
diferentes entre las metodologías. SoilPar 


0.35 A 
| | O | 
030 q o o o 
=— Rosetta : : : : ñ : 
1 E Rosetta SSCB 
) Rosetta SSC 
0.25 === 


A 
A  SoilPar 
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0.20 


Parámetro 0, (adim) 
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Figura 2. Comparación de la estimación de la humedad residual Or mediante Rosetta SSCB 88, 50 


SoilPar y ecuaciones del cuadro 1. 
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Figura 3. Comparación de la estimación de la humedad saturante Os mediante Rosetta SSCB 0,,-8, ¿y 
SoilPar y ecuaciones del cuadro 1. 
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Figura 5. Comparación de la estimación del parámetro n mediante Rosetta SSCB 0,,-8, ..y 
SoilPar y ecuaciones del cuadro 1. 
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Figura 6. Comparación de la estimación del parámetro Ks mediante Rosetta SSCB 8,,-0, ¿y 
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tiende a valores menores que Rosetta, mientras 
que las ecuaciones de regresión (Zimmermann 
y Basile 2007) tienden a valores mayores en un 
orden de magnitud. Para aplicar estas últimas 
a las texturas arcilla limosa y franco limosa 
(cuadro 1), se asumió que la densidad aparente 
DA es igual a la densidad húmeda Ó,, que la 
porosidad y es igual a 6, y que la densidad seca 
S, es igual a 6, / 1.25 (es decir, se consideró una 
humedad del 25% promedio en peso). Para 
la textura franco arcillo limosa no fue posible 
aplicar las ecuaciones debido a la insuficiente 
información granulométrica. 

En la figura 5 se observa una tendencia 
inversa para la estimación del parámetro n, 
es decir, el software SoilPar tiende a valores 
mayores de n respecto a Rosetta y las demás 
estimaciones lo hacen por defecto. En este caso, 
las diferencias entre los distintos métodos son 
de menor magnitud que en el caso anterior. 

Con respecto al pronóstico de la con- 
ductividad hidráulica saturada (figura 6), 
las tendencias son muy disímiles entre los 
métodos evaluados. Incluso entre las ecua- 
ciones de Puckett y Jabro incluidas en el 
software SoilPar, llegan a valores de un orden de 
magnitud diferente para una misma muestra 
de suelos. 

Esto da la pauta de que se debe considerar 
con precaución el empleo de estos programas, 
ya que las ecuaciones ajustadas internamente 
dan resultados muy distintos para una misma 
caracterización físico-química de un suelo. 


Aplicación de FPTs en la cuenca del arroyo 
Ludueña 


Ante esta diversidad de resultados, es necesa- 
rio contar con elementos dejuicio que permitan 

FPTs 
comparativos 


seleccionar las más convenientes. 


Los valores disponibles 
consistieron en parámetros calibrados con 
el modelo hidrológico-hidráulico SHALL3 
(Zimmermann y Riccardi, 2000). En el modelo 
se contemplan las interacciones entre los 
procesos de superficie en la zona no saturada 


(ZNS) y en la zona saturada. Los flujos en 


la ZNS se simulan mediante la ecuación de 
Richards, contemplando el modelo de Brooks 
y Corey para la curva de retención. El modelo 
SHALLS3 fue puesto en operación en el sistema 
hidrológico del arroyo Ludueña (Santa Fe, 
Argentina), y posteriormente fue calibrado 
utilizando información freatimétrica, pluvio- 
gráfica y limnigráfica en el periodo 1994- 
1996 (Zimmermann, 2003). Dado que los 
parámetros se ajustaron en forma global 
para la cuenca, sólo constituyen valores de 
referencia generales. 

Las FTPs analizadas se han aplicado a las 
series de suelos presentes en la cuenca del A? 
Ludueña, cuyas asociaciones se combinan 
arealmente de acuerdo con la distribución que 
se muestra gráficamente en la figura 7. 

Los parámetros se han promediado por 
horizonte (A y B, según cartas de suelo de INTA) 
y ponderado arealmente de acuerdo con los 
porcentajes de asociaciones y series analizadas. 
La metodología de ponderación para obtener 
parámetros por serie de suelos, partiendo de 
los parámetros calculados para cada horizonte, 
se basó en los siguientes criterios: (a) para la 
conductividad hidráulica vertical se calculó 
la media armónica de las conductividades de 
cada estrato (ecuación (4)) y (b) para el resto de 
los parámetros, el promedio ponderado por el 
espesor de cada estrato (ecuación (5)): 


K= (4) 
e: 
=1 


1 


i 


1 


¿Y 
X= 9s se (5) 


donde X,, es el valor medio de un parámetro 
cualquiera; x, el valor del parámetro para el 
espesor e; N, el total de horizontes de la serie; 
K,, la conductividad hidráulica media; K, 


conductividad del horizonte, y L es el espesor 
N 
total de la serie L= ys Es 


i=1 
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Figura 7. Distribución de las asociaciones de suelos en la 
cuenca del arroyo Ludueña, Santa Fe, Argentina. 


Las series puras de suelos mencionadas 
precedentemente se combinan entre sí, en 
diferentes proporciones, dando lugar a aso- 
ciaciones de series de suelos. Para estimar los 
parámetros de las asociaciones se adoptó una 
metodología de ponderación basada en el 
porcentaje de ocupación areal de cada serie 
dentro de cada asociación y las propiedades 
correspondientes a cada serie: 


M 
Y x, A; (6) 


j=1 


Xas= 


>| 


donde X, es el valor medio de un parámetro 
cualquiera para la asociación; Xy el valor del 
parámetro para la serie de suelos j; M, el total 
de series que componen la asociación, A, el área 
de ocupación total de la asociación, y A, es el 
área de ocupación de la serie ¡ en la asociación 
en cuestión. 

El cuadro 3 presenta los valores ponderados 
arealmente y promediados por horizontes de 
los parámetros analizados para las diferentes 
En el 
cuadro se han resaltado las formulaciones que 


FPTs consideradas en los análisis. 


aproximan en +10% los valores calibrados. 


Puede observarse que en cuanto a la 
conductividad hidráulica saturada, los valores 
propuestos por Saxton et al. (1986) y Rosetta 
con la información completa están dentro del 
entorno propuesto. 

En cuanto al parámetro az, la formulación de 
Wosten et al. (1999) y Rosetta con la información 
completa también se encuentran dentro del 
entorno propuesto. 

Para pronóstico del parámetro n, el progra- 
ma Rosetta, en todas las versiones, da mejores 
resultados que el resto de las formulaciones. 

El parámetro que presenta mayor sensi- 
bilidad al proceso de infiltración en el 
modelo hidrológico es el K,, mientras que los 
parámetros que regulan la curva de retención 
de humedad tienen un rol secundario debido 
a que la variación de humedades durante las 
simulaciones no es muy grande en la zona de 
estudio considerada. 

Los resultados alcanzados son promisorios 
y validan la aplicación conjunta de la 
metodología de ponderación propuesta y el 
software Rosetta para una buena estimación 
de parámetros hidráulicos. Eventualmente 
pueden utilizarse las formulaciones de Saxton 
et al. (1986) y Wosten et al. (1999) para estimar 
la conductividad hidráulica saturada y el 
parámetro Q. de la ecuación de van Genuchten, 
respectivamente. 


Conclusiones 


Se han utilizado diferentes procedimientos 
para estimar propiedades hidráulicas de 
suelos limosos, representativos de un sector 
de la partiendo de 
información granulométrica y de parámetros 


llanura argentina, 
físico-químicos disponibles en las cartas de 
suelos de INTA. Se emplean funciones de 
pedotransferencia (FPT) incluidas en los 
códigos SoilPar y Rosetta, junto con ecuaciones 
de regresión entre parámetros hidráulicos y 
parámetros físicos dadas por un conjunto de 
autores. 

Para el conjunto de las muestras analizadas, 
el análisis individual de cada una reveló 
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Cuadro 3. Valores comparativos de parámetros de van Genuchten y conductividades hidráulicas saturadas 
obtenidos para las texturas analizadas en la zona de estudio. 


K 
Autores tenh) 
Horizonte “A” Horizonte “B” Promedio armónico 

SoilPar (ec. Jabro) 0.836 1.141 0.965 
SoilPar (ec. Puckett) 2.141 0.011 0.021 
Rosetta SSCB 6,.-8, ., (cuenca Ludueña) 0.355 0.190 0.248 
Rosetta SSCB (cuenca Ludueña) 0.474 0.376 0.420 
Rosetta SSC (cuenca Ludueña) 0.498 0.714 0.586 
Vereecken et al., 1990 3.587 6.200 4,545 
Cosby et al., 1984 0.541 0.334 0.413 
Wosten et al., 1999 3.213 2.785 2.984 
Saxton et al., 1986 0.428 0.063 0.109 
Rawls y Brak, 1985 0.066 0.011 0.019 
Campbell y Shiozawa, 1994 0.013 0.000 - 

Jabro, 1992 0.917 1.251 1.058 
Jaynes y Tayler, 1984 0.616 0.388 0.476 
Gulser y Candemir, 2008 5.668 1.876 2.819 
Wosten et al., 2001 153 1867 1.446 
Puckett, 1995 2.142 0.011 0.021 


Parámetros calibrados en la cuenca del A* Ludueña 


(Zimmermann 2003) 0.125 (0.225-0.030) 


Autores a lea 
Horizonte “A” Horizonte “B” Promedio 

Carsel y Parrish, 1988 0.0200 0.0050 0.0125 
SoilPar 0.0015 0.0010 0.0013 
Rosetta SSCB 6,,-8, .,,, (cuenca Ludueña) 0.0030 0.0039 0.0035 
Rosetta SSCB (cuenca Ludueña) 0.0063 0.0141 0.0102 
Rosetta SSC (cuenca Ludueña) 0.0066 0.0142 0.0104 
Wosten et al., 1999 0.0076 0.0036 0.0056 

Rawls y Brak, 1985 0.0167 0.0140 0.0154 
Zimmermann, 2006b (IIS) 0.0293 0.0763 
Parámetros calibrados en la cuenca del A* Ludueña 0.0040 

(Zimmermann, 2003) 

n 
Autores G 
Horizonte “A” Horizonte “B” Promedio 

Carsel y Parrish, 1988 1.41 1.09 1125 

SoilPar 1.68 1.14 1.41 

Rosetta SSCB 0,.-6,., (cuenca Ludueña) 1.56 152 1.44 

Rosetta SSCB (cuenca Ludueña) 1.60 1.34 1.47 

Rosetta SSC (cuenca Ludueña) 1.59 1.34 1.47 
Wosten et al., 1999 1.20 1.09 1:15 

Rawls y Brak, 1985 1.30 1.20 1.25 
Zimmermann, 2006b 1.38 1.26 1,32 
Parámetros calibrados en la cuenca del A* Ludueña 

(Zimmermann, 2003) 1.56 
una gran dispersión entre los diferentes precaución, analizando las compatibilidades 
procedimientos empleados. Esto da la pauta entre las muestras de suelo de estudio y las 
de que los predictores se deben emplear con que se utilizaron originalmente en la propuesta 
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de las ecuaciones intrínsecas de cada FPT. Es 
decir, hasta el momento no pueden obtenerse 
ecuaciones aplicables “universalmente”. 

Se estimaron los parámetros medios 
ponderados para la cuenca del A” Ludueña 
y los pronósticos se compararon con valores 
obtenidos mediante calibración de un modelo 
matemático hidrológico (Zimmermann, 2003). 

A partir de los resultados obtenidos en 
este trabajo pudo concluirse que entre los 
códigos empleados, el que mejor se aproximó 
a los parámetros hidráulicos calibrados para 
la zona de estudio fue el programa Rosetta, 
suministrándole la información completa 
(composición textural, densidad aparente y 
dos puntos de la curva de retención: punto 
de marchitez permanente y capacidad de 
campo). Eventualmente pueden utilizarse las 
formulaciones de Saxton et al. (1986) y Wosten 
et al. (1999), para estimar la conductividad 
hidráulica saturada y el parámetro a. de la 
ecuación de van Genuchten, respectivamente. 

Los resultados alcanzados al presente son 
promisorios y validan esta metodología de 
ponderación sencilla para ser aplicada en otras 
cuencas de la región. 


Glosario de variables 


A (%): 

C (%): contenido de arcilla. 

CIC (cmol (+) kg”): capacidad de intercambio 
catiónico. 


contenido de arena. 


CO (%): contenido de carbono orgánico. 
ev: coeficiente de variación geométrico. 
d: diámetro geométrico. 

DA (gcm”): densidad aparente. 

K (cm h*): conductividad hidráulica saturada. 
L(%): contenido de limo. 

MO (%) contenido de materia orgánica. 


m (adimensional): parámetro empírico del 
modelo de Van Genuchten. 
n (adimensional): parámetro empírico del 
modelo de Van Genuchten. 
concentración de hidrogeniones. 
variable cualitativa binaria (0, 1), 
que indica si el estrato es superficial 
(1) o no (0). 


pH: 
top: 


al (cm): parámetro empírico del modelo de 
Van Genuchten. 
P.: curtosis adimensional. 


S, (g/cm*): densidad húmeda. 
9, (g/cm?): densidad seca. 

n (adimensional): porosidad. 
y (kPa): 


m 


succión correspondiente a un con- 

tenido de humedad determinado. 

9, (m* m>”): contenido de humedad a -20 kPa. 

9 ,, (m* m”): contenido de humedad a -33 kPa, 
capacidad de campo. 

9, (m'* m”): contenido de humedad a -40 kPa. 

9 y, (mé m?): contenido de humedad a -200 kPa. 


9, ¿y (m* m”): contenido de humedad a -1 
500 kPa, punto de marchitez 
permanente. 

9, (m* m?): contenido de humedad residual. 

9, (m*m”): contenido de humedad para el 


suelo saturado. 


Los parámetros granulométricos se obtie- 
nen de la siguiente forma: 

Escala sedimentológica b: 4 = log, d, siendo 
d el diámetro de la partícula. Los parámetros 
estadísticos considerados para las distribu- 
ciones granulométricas son diámetro medio 
0,, desvío estándar ou, curtosis K, curtosis 
adimensional f,. Los mismos se calculan en la 
escala ( de la siguiente manera: 


N N 1/2 
0 = Y if, o= S (5-0) | 
i=l 3 
N di d 
ez (o, = Om Yi B, = o! , donde f, es la 


fracción en peso asociada con el diámetro d.. 
El diámetro medio geométrico d,, el desvío 
estándar geométrico o, y el coeficiente de 
variación geométrico CV, se obtienen a través 
de los respectivos valores definidos en la 
escala q mediante las ecuaciones: 


d 
d¿=2 0,=2" CV, = MEN 
, , (0) 
ES 


Recibido: 06/08/09 
Aprobado: 10/08/10 
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Abstract 


ZIMMERMANN, E.D. és BASILE, P.A. Estimation of hydraulic parameters in silty soils 
using different pedotransfer functions. Water Technology and Sciences, formerly Hydraulic 
engineering in Mexico (in Spanish). Vol. 1, No. 1, January-March, 2011, pp. 99-116. 


This paper presents different procedures that have been used to estimate hydraulic properties of 
silt soils, representative ofa sector ofthe Argentinean flatlands. Grain size, and physicochemical 
information of soil maps printed by the National Agriculture Technology Institute (INTA) were 
used. Pedotransfer functions (PTE) included in the SOILPAR codes (Acutis and Donatelli, 
2003) and Rosetta (Schaap et al., 1999) together with empirical relationships proposed by 
other authors were used. For the set of 136 analyzed samples of soils, the application of the 
procedures showed an important dispersion of results. Consequently, the PTF applied must 
be used with caution, and the compatibilities between the studied soil samples and those that 
were used originally in the proposed equations must be analyzed. The average parameters 
estimated in the Ludueña basin (Santa Fe, Argentina) were compared with values obtained 
by calibration of a hydrologic mathematical model. The results showed that the Rosetta code 
with complete information (textural composition, bulk density and two points of the retention 
curve) had better prediction of hydraulic parameters than the other methods. The formulations 
of Saxton et al. (1986) for saturated hydraulic conductivity and Wosten et al. (1999) for van 
Genuchten's parameter Q. can be alternatively used to the Rosetta model. The results reached 
to date are promissory, and validate the methodology for its application in other regional river 
basins. 


Keywords: silty soils, hydraulic parameter, pedotransfer functions, flatlands of Argentina. 
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